Preprint typeset in JHEP style - HYPER VERSION 



Fermion Families from Two Layer Warped Extra 
Dimensions 



Zhi-qiang Guo^ and Bo-Qiang Ma^* 

^ School of Physics and State Key Laboratory of Nuclear Physics and Technology, Peking 
University, Beijing 100871, China 



Abstract: In extra dimensions, the quark and lepton mass hierarchy can be reproduced from 
the same order bulk mass parameters, and standard model fermion families can be generated 
from one generation in the high dimensional space. We try to explain the origin of the same 
order bulk mass parameters and address the family replication puzzle simultaneously. We 
show that they correlate with each other. We construct models that families are generated 
from extra dimensional space, and in the meantime the bulk mass parameters of same order 
emerge naturally. The interesting point is that the bulk mass parameters, which are in same 
order, correspond to the eigenvalues of a Schrodinger-like equation. We also discuss the problem 
existing in this approach. 



Keywords: Fermion Families, Warped Extra Dimensions, Mass Hierarchy 



*Email:mabq@phy. pku.edu.cn 



Contents 

|l]. Introduction [l] 

|2]. The Setup || 

|3|. Examples of metric for finite generations 

^ Further Discussions |l3| 

[H]. Conclusions 15 

\A[ Explicit solutions for massive modes and zero mode: metric I 15 

1^. Numerical examples for finite generations 17 

[^ . Property of hypergeometrical function 18 

|D| . Solutions for zero mode and massive modes: metric II 18 

|E| . Another metric example for orthogonality 19 

Ricci scalar curvature for metric 20 



1. Introduction 

In recent years, extra space dimensions have played an important role in particle physics, 
gravity and cosmology. Many old problems get elegant solutions in extra space dimensional 
background with new perspectives [^, ^. The extra space dimensions have been used to 
address the problems in flavor physics: why the masses of quarks and leptons distribute in 
a large range and have obvious hierarchy structure (the fermion mass hierarchy puzzle); why 
the heavier generations replicate the lightest generation with the almost same properties (the 
fermion family replication puzzle). There have been very interesting progress that both of 
them get answers in several different approaches ^, ^, ^, • 

In extra dimension, the fermion mass hierarchy has geometrical origin. It arises from the 
small overlap of wave functions in the extra space dimensions. The hierarchy mass structure 
can be reproduced with the bulk mass parameters of the same order in 5-dimension warped 
space (Randall-Sundrum model) successfully (for several numerical examples, see |Q).^ How- 
ever, a new puzzle arises naturally: why the 5-dimension bulk mass parameters are of the same 

^For definiteness, our following discussions will be based on the concise numerical example given in and 
we hope that our discussions can apply to other examples by some modification. 
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order, or whether we can give an explanation to the origin of the same order 5-dimension bulk 
mass parameters. It is easy to understand that this new puzzle correlates closely with the 
fermion family replication problem. One bulk mass parameter stands for one fermion family 
or one fermion flavor, so to explain the origin of the same order bulk mass parameter is doing 
the same thing to explain the family replication problem. 

This new puzzle can be addressed in the approach that the three generations of standard 
model (SM) can be generated from one generation in the high dimensional space. This approach 
has been adopted in several papers in which the 6-dimension spacetime is reduced to 

4-dimension spacetime directly. It is found that the SM families correspond to the zero modes 
of the high dimensional equation of motion. In these approaches, one missed the chance to 
give answers for the origin of the same order bulk mass parameters in 5-dimensions. In the 
present paper, we will adopt an alternative approach. We suggest a 6-dimension metric Ansatz 
of special two layer warped structure, 

ds"^ = B{zf [A{yf {rif.^dxf'dx'' + dy^) + dz'^] . (1.1) 

As we will show below, with the help of the special structure of this two layer warped metric, 
we can reduce the 6-dimension (4-l-l)-|-l spacetime to 5-dimension 4-1-1 spacetime at the first 
step. We found that the induced 5-dimension equation of motion for fermions is similar to 
that analyzed in IH, while the bulk mass parameter in this induced 5-dimension equation 
corresponds to the eigenvalue of a 1-dimension Schrodinger-like equation, which is conformed by 
the fermion wave function in the sixth space dimension. Yet, because the induced 5-dimension 
equation of motion is similar to that analyzed in [^, we might expect that this induced 5- 
dimension equation can be applied in this model by some modification. Hence we reduce the 
problem of the same order 5-dimension bulk mass parameter to the eigenvalue problem of 
a second order differential equation. Because the eigenvalues of a 1-dimension Schrodinger- 
like equation are of the same order generally, we see that the same order 5-dimension mass 
parameters emerge naturally. If we further reduce the 5-dimension spacetime to 4-dimension 
spacetime, with the help of the model [^, the same order 5-dimension bulk mass parameters 
will produce the hierarchy structure of the 4-dimension physical fermion mass. By this new 
approach, we can address the new puzzle suggested above. We give an explanation for the 
origin of the same order bulk mass parameters. Because one 5-dimension bulk mass parameter 
corresponds to one family, we provide an answer for the origin of the families simultaneously. 
However, there exist some problems in this approach, we will discuss these problems in detail 
and suggest methods to bypass them. 

According to the approach suggested above, we construct a model that the 5-dimension 
bulk mass parameters of the same order emerge naturally, and hence the standard model 
families are generated simultaneously. In fact, we can show that the 5-dimension bulk mass 
parameters correspond to the eigenvalues of a second order ordinary differential equation of 
the Sturm-Liouville type (similar to the Schrodinger equation in 1-dimension). Therefore, 
like the energy spectrum of Hydrogen atom, the bulk mass parameters or the families can be 
generated from the eigenvalues of Schrodinger-like equation with proper boundary conditions. 
The details will be introduced in Sec. 0. Several examples are given in Sec. I We give further 
discussions and conclusions in Sec. ^ and Sec. ^. 
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2. The Setup 



In this section, we introduce the setup in detail. We start with the action of a bulk Dirac 
fermion in six dimension spacetime. The metric Ansatz for spacetime takes the form 

ds^ = B{zf [A{yf {rjf.^dxf'dx'' + dy"^) + dz"^] , (2.1) 

where r/^,y = diag(— 1, 1, 1, 1). Note that the metric has the special two layer warped structure. 
By this special metric structure, one can reduce the 6-dimension spacetime to 5-dimension 
spacetime at the first step, then further reduce the 5-dimension spacetime to 4-dimension 
physical spacetime, as we will introduce in detail below. We suppose that the extra dimensions 
both are intervals.^ 

The bulk action for this fermion is given by the usual form, 

S = j d^xdydz^ c^^VDm^ - Dm^ ei^V^i'] -im^^^, (2.2) 

where e^^ is the sechsbien, and Dm = Qm + \^fi^ah, ^ah = \ \^aj ^h] is the covariant derivative 
of spinor in curved spacetime. a and M = 0,1,2,3,5,6 stand for the flat spacetime indices 
and the curved spacetime indices respectively. The Dirac equation in 6-dimension spacetime 
requires m to be a real number. We choose the gamma matrices representation as follows. 




where /i = 0, 1, 2, 3 and fi* are the usual Pauli matrices. 
By use of the metric Ansatz, the action reduces to 



S = Si + S5 + Se- j d'^xdydz^^i m^^, (2.4) 
S4 = y d'^xdydz^i^^Bizy^Aiyy^ [^T^'d^^ - d^^T^^ 
S^ = j d^xdydz^ I '-B{z)~^A{y)-^ r^^s^ - d^^'T^'^] 
Se = j d^xdydz^/^ !^^B{zy^ [^T^de^^ - Se'lT^^] | . 
Varying the action with respect to ^, we obtain the equation of motion, 

A(y)-i + (95 + 2A-^A') ^} + r*^ (^de + ^ -mB^ = 0, (2.5) 

^The conventional way in extra dimensions is to suppose the extra dimensions as orbifolds, and that the 
boundary conditions are determined to be the Israel's junction conditions as adopted in j^. However, 
as suggested in the interval approach is more convenient in some cases, and it can produce more general 
boundary conditions. In our work, we found that it is necessary to adopt the interval approach, at least when 
we deal with the boundary conditions in the sixth space dimension. Of course, it is also more convenient. More 
about the interval approach, see [Q. 
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where A' = '^^}y^ , B = , = dy,dQ = dz, and the boundary term, 

55'bound = '^'55 bound bound) 

(2.6) 

555bound = -lj <^^^dz B-^A-H^T''^''^ , 

556bound = -'2 J d^^dy S-i^^r^^jf , 

where we denote by [X]^ the quantity — X\q. Denoting ^' = [ ) , where xi and X2 are 
four-component Dirac spinors, we rewrite Eq. ( p. 51) as 

A{y)-' {j^^df,X2 + l' {85 + 2A~'A') X2} + (^e + ^B-^B^ xi-mBxi = 0, (2.7) 

A{y)-' {j^'d^Xi + lHd5 + 2A-'A')xi}+ -(^Oe + ^B-^B^ X2-mBx2 = 0- (2.8) 

Now we make the conventional Kluza-Klein (KK) decomposition. We expand xi X2 
with spinor 'ilj{x'^,y) in 5-dimension spacetime as 

Xiix'', y,z) = Y,Fn{z)^n{x^,y), X2(x^,y,z) = ^ G„(z) V^„(x^, y), (2.9) 

n n 

in which ipn{x^^y) conforms to the Dirac equation in 5-dimension spacetime, 

A{y)-^ {^d^Mx^. y) + 7' (^5 + 271-1.4') y)} - A,Vn(x^ v) = 0- (2.10) 

As in 6-dimension spacetime, the Dirac equation in 5-dimension spacetime requires A„ to be 
real numbers. This can be verified by multiplying the two sides of Eq. ( 2.1C| ) by ipnix'^,y)- 



Here we note that it is critical that A„ must be real numbers, as it will be obvious in the 

following discussions. With the help of Eq. ( p.lO| ), Eqs. (|2^)-(|2^ can be solved by the 
following Ansatz, 

^ + ^B-^B^ Fn{z) - mBFn{z) + A„G„(z) = 0, (2.11) 

{j'z ^ ^"^^^ + rnBGniz) - A„F„(z) = 0. (2.12) 
These equations can be simplified further by the transformations, 

Fn{z) = B{z)-'Fn{z), Gniz) = B{z)-'Gn{z), (2.13) 
in which e = |- Then we obtain the equations 

^ - mB^ Fn{z) + XnGn{z) = 0, (2.14) 

l- + mB\Gn{z)-XnFn{z) = Q. (2.15) 
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For a zero mode (A = 0), these bulk equations decouple and are easy to be solved. The 
solutions are given by 



Fo{z) = ^=e^p ( [ mB{C)dc] or 0,Go{z) = ^=exp ( - [ mB{C)dc] or 0.(2.16) 
V'^VO \Jzo J JlNr. ^ ■''^0 J 



We introduce / of the length dimension in order to make the normalization constants to be 
dimensionless. It will become explicit in examples in the next section. For the massive modes, 
we can combine the first order differential equation to obtain second order equations 



dz^ 



Fn{z) + 
Gn{z) + 



mB — m?B^ 
mB — m^B^ 



Fn{z) + XlFn{z) = 0, 
Gn{z) + XlGniz) = 0. 



(2.17) 
(2.18) 



Rewriting them in another form, we see that they are similar to the one dimensional Schrodinger 
equations 



- ^Fn{z) + V{z)Fn{z) = XlFniz), 

- ^Gn{z) + V{z)Gn{z) = XlGniz), 



with potentials 



V{z) =mB + ra^B"^, V{z) = -mB + m^B 



2d2 



(2.19) 
(2.20) 

(2.21) 



In the following discussions we will illuminate that such a setup gives answers to the puzzle 
we proposed in the introduction. According to the setup above, we realize an interesting fact 
that Eq. (2.10) is similar to the equation analyzed in the model if we choose A{y) to be a 
slice of the anti-de Sitter (AdS) metric, i.e., the RS spacetime. The differences are that there is 
gauge field background in the model , and that in that work the extra dimension is adopted 
to be an orbifold. However, for fermions, the equations of motion in these two cases are of 
almost similar features, and the gauge field background only makes the boundary conditions 
more involved. It is not difficult to add the gauge field background as in the model to the 
above setup. When we further reduce the 5-dimension spacetime to 4-dimension spacetime, we 
can advance with the help of the model |^]. A„, in the above setup corresponds to the bulk mass 
parameter in 5-dimension spacetime in the model |^]. Of course, it is obvious that one bulk 
mass parameter corresponds to one family in 5-dimensions. Now we can understand how the 
same order bulk mass parameters or families are generated from extra space dimensions. The 
bulk mass parameters are eigenvalues of Schrodinger-like equations ( |2.19D -( p.20 ), so generally 
they should be of the same order. The eigenstates of equations ( |2.1£ )-( [2.20 ), belonging to 
the eigenvalues A„, correspond to the generations in 5-dimensions. When we reduce further 
the 5-dimension spacetime to the physical 4-dimensions, these generations in 5-dimensions can 
produce the generations in 4-dimension spacetime. So the families in physical 4-dimension 
spacetime are generated simultaneously. Of course, the eigenstates should be normalizable in 
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order that we can get the effective 5-dimension action after integrating out the sixth dimen- 
sion. We will discuss the normalization conditions in the next section. However, a problem 
arises immediately from the following contradiction: On one side, the eigenvalue problem of 
Eqs. ( p.l9 )-( 2.20| ) is of the Sturm-Liouville type. The characters of the Sturm-Liouville eigen- 
value problem are that the number of eigenvalues is infinite and the size of eigenvalues is 
non-bounded, i.e., the eigenvalue series becomes large monotonously; on the other side, as it 
has been illuminated obviously in the papers [^, the larger bulk mass parameters produce 
lighter fermions mass in 4-dimensions. Therefore, the eigenstates of Eqs. ( 2.19D -( p.20 ) produce 
infinite light fermion generations. However, no lighter generations are discovered by experi- 
ments so far. Hence we need another mechanics to cut off the infinite series and select only 
several eigenstates. The left eigenstates correspond to generations in 4-dimensions. 

There also exists a problem about the zero mode ( p. 16 ). By the numerical examples in Q, 
the zero bulk mass parameter in 5-dimensions produces a very heavy fermion in 4-dimensions, 
and it is heavier than the SM generations. So it does not correspond to the physical generations. 
If the zero mode is permitted by the boundary conditions and the normalization conditions in 
our model, there would exist a generation that has not been discovered by experiment so far. 
However, the zero mode is more subtle in the example we will discussed in the next section. 
We will discuss this problem in more detail in that section. 

Now we suggest several approaches to deal with the problem about the infinite eigenvalues. 

(1) An immediate proposal is that one chooses a 6-dimension spacetime in which the sixth 
dimension is not continuous but discrete. For example, if we discrete the finite interval to be 
finite points, the induced Eqs. ( 2.19D -( p.20 ) will be difference equations. The number of their 



eigenvalues is finite naturally. There has been a similar investigation for gravity, see |13|. 

(2) Another bizarre proposal is that we suppose the sixth dimension to be timelike. In the 
Ansatz Eq. ( |2.1| ), we have chosen the sixth dimension to be spacelike. Instead we can choose 
it to be timelike. This leads to the metric with two time dimensions [|l4|.^ The alternative 
metric Ansatz is 

(2.22) 



ds"^ = B{zf [A{yf ['n^^dx^'dx" + dy'^) - dz' 



In Eq. (|2j), we let 
produces the equations. 



, with others keeping invariant. The same procedure 



-XriB ^B — A? 



;Gn{z) - B-'B^Gn{z) + m^Gn{z) + 



XnB ^B — Xi 



dz'^ ' dz 
We give a simple example in which B{z) = constant. The solutions are 

F{z) = Ci exp'''' + C2 exp-^'=^ k = Vm^ - X^. 



Fn{z) = 0, 
GJz) = {). 



(2.23) 
(2.24) 

(2.25) 



Here we have omitted the subscripts. If we impose the boundary conditions, 

F|o = 0, F\r = {), 



'For extensive investigations on two-time physics, see ||l5 



(2.26) 



- 6 - 



then A must conform to 



kR = - X'^R = nvr, n = 1, 2, 3, • • • . (2.27) 



As we emphasized above, A and m must both be real numbers. Eq. ( p. 27 ) has solutions only 



for finite natural number. The number of the eigenvalues depends on the size of m, hence there 
are a finite number of eigenvalues. 

(3) However, in the following section we will adopt another approach, i.e., we can obtain a 
finite number of eigenstates by choosing the metric B{z) delicately. Although B{z) produces 
infinite eigenstates generally, it is possible that some of them can produce only finite eigenstates. 
We will focus on this possibility in the following section, and give concrete examples for finite 
generations. 



3. Examples of metric for finite generations 



In this section, we will suggest metric examples which can produce finite generations. We 
discuss the appropriate normalization conditions for Fn{z), Gn{z) at first, then we analyze 
examples in detail. 

We begin with the action (U). With the help of Eqs. (|2.10D , and ([2l^ ), the action 

(|2.4|) can be rewritten as 



S = d xdy K„ 



dzB^ {F:^Fn + Gl,Gn 



dzB'' 



TP* TP J- r^* c 



dz {F^Ffi -\- GyjiGfi) , 



\ A„ + AJ^ 




/ 2 





An + A^ 



(3.1) 
(3.2) 
.(3.3) 



Notice that ^' = ^-tpO = (x2,Xi)- Eqs. (|1[]3]|) are satisfied for all modes, including zero 
modes and all massive modes. 

In order to get the conventional effective 5-dimensional action 

^5eff = ^ j d^xdy i^^A^ [^nT^^sV'n - d^i)nl^i^n + tpnl'^d^tpn " S^V'nT'^V'n] | 

- ^ / d^xdyA^iXnMn, (3.4) 

we consider two cases below: 

Case (I): The first case is that the normalization conditions 



Kmn = Jdz (F^Fn + Gl.Gn) = 5„ 



(3.5) 
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are satisfied. As in the standard Sturm-Liouville case, we can convert the normahzation 
conditions ( |3.5| ) to the boundary conditions. By Eqs. ( p.l7| ) and (2.18), we have 

f-R' 



JR 



F —F* 



F* —F 



G^-Gr 



\R' 
\r 



(A„ + X^){F^Gr. 



'^m^nj \r ■ 



(3.6) 
(3.7) 



In the last line, we have used the bulk equations (2.14) and (2.15) to simphfy these expressions. 
We can also get ( |3.7| ) from Eqs. ( |2.14| ) and (2.15) directly. There are two types of concise choices 
to make the normahzation conditions satisfied, 



{a): F\r 
or (6) : F\r 



0, F\r, = 0; 

G\r, F\ri = G\ri 



(3.8) 
(3.9) 



Then for real Am, the orthogonahty is ensured by appropriate boundary conditions. In this 
case, we can get Eq. ( |3.4[ ) from Eq. ( |3.1| ) via Eq. ( |3.5| ) in a straight way. 

Case (II): The second case is that the normahzation conditions Eq. ( p.5p are not satisfied. 
In this case, K and M are both matrices, which means that different KK modes are mixed 
not just among the mass terms, but also among the kinetic terms. At the first sight, it seems 
that we can not get the the conventional effective 5-dimensional action Eq. ( |3.4| ). However, if 
K is positive-definite and the number of KK modes is finite^, we can redefine the fermion field 
to get an action, which has the same form with that of Eq. (|3.4| ). The difference is that the 
eigenvalues A„ are modified to different size. From Eq. ( [3.2| ) and Eq. ( |3.3D , we know that K 
and M are both hermitian. A positive-definite hermitian matrix K can be diagonalized as 



K = yUy = h'^h, h = VXv, 

A = diag(Ai,A2,---,A„), 
VX = diag(v/Ar, ^/A^, y/K). 



(3.10) 



In the above expressions, Aj > 0, i = 1, 2, • • • , n, as we have supposed that K is positive- 
definite. Redefine t/'n as 



then in the new basis ■(/'n, M becomes 

After diagonalizing M by ?7, we have 

M = U'^AU, 

A = diag(Ai, A2, • • • , A„). 



(3.11) 



(3.12) 



(3.13) 



*We will give numerical examples to show that such conditions can be satisfied in Appendix I 
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The action ( ^.l] ) can be reduced to the form like that of action (| 



(3.14) 



In the above, we have given the normahzation conditions. The criteria are that we can 
integrate out the sixth dimension to get an effective 5-dimension action. Now we suggest an 
example that can produce finite generations. In Eqs. ( 2.17] ) and ( |2.18| ), we suppose that 



B(z) 



+ a 



+ b 



s, a, 6, a; > 0. 



, U, Hv, v^, 



(3.15) 



As in the models lo can be regarded as the characteristic energy scale of the sixth dimension. 
We will see that it determines the size of KK modes below. The role of the dimensionless 
parameters s and b will become obvious after we give the solutions of Eqs. ( p.l7| ) and ( |2.18| ). 
The conditions o, 6 > ensure that the metric is well behaved in the interval (— oo, oo). 
Eqs. (2.1?) and ( p.l8| ) can be solved by hyper geometrical functions, 



F{z) = Cle-^"^''{e'^'' + 6)^~^hypergeom [ p - fi + - p - n + v;! - 2^ 



+ C2e'^'^^(e'^^ + 6)-^-''hypergeom [ p + p + u,l - p + p + u-1 + 2p 



oU)Z 



+ h 



,(3.16) 



where p = ^s{l - I), p = ^ [^sY (f)^ - (A)^ and u = ^[{^' - (A)^ C, and C, are 
constants. For the sake of simplicity, we omit the subscript n. We only display the solution 
for F{z) explicitly. The solution for G{z) can be determined by F{z) through Eq. ( 2.14] ) 
or by Eq. ( 2.1^ ) directly. Now let us investigate this solution. For hyper geometrical func- 
tion^ Hypergeom(a, /3; 7, ^), when Re(7 — a — (3) < 0, it diverges at = 1. In the solution 
(3.16), we have 



Re(7 -a- (5) = Re{-2u) < 0, 



(3.17) 



so when z 



oo, i 



'-+b 



1, a singularity happens. If we choose z to be a finite 



interval [i?, i?'], then the solution ( 3.16 ) is well behaved in this range. Imposing the boundary 
conditions ( ]3.8] ) or ( |3.9[ ), we get infinite eigenvalues generally. However, we find that the 
following choice can produce a finite number of eigenvalues. Given that z to be a semi- 
infinite interval oo), then the solution ( 3.16] ) develops a singularity when z ^ oo. This 
singularity makes the integral in ( |3.5D to be divergent. In order to make the integral to be 
finite, the hypergeometrical series ( ]3.16] ) must be cut off to be a polynomial by the requirement 
a = —{n — 1), or (3 = —{n — 1), n = 1, 2, 3, • • •. In Eq. ( 3.16) ), we choose 



1 — p — p + u = —{n — 1), n = l,2,3,' 



(3.18) 



^ About the property of hypergeometrical function, see Appendix |^. 



- 9 - 



Because A and m are real numbers, as we have discussed above, Eq. ( 3.18 ) might have solutions 
if 



1 < n < 



m 

— s. 

UJ 



(3.19) 



Obviously, the size of n are limited by the parameters in the metric, then only finite eigenvalues 
are permitted. We see that it is important that A and m are real numbers again. The conditions 
(3.18) are required by the boundary conditions when z — > oo. This boundary condition restricts 
the solutions to the form ( |A.l| )-( |Ar3| ). We display these solutions in Appendix ^ explicitly. We 
see that these solutions are determined completely up to normalization constants. For these 
solutions, Fn{z), Gn{z) x~'^" in the symbols in Appendix^, when z — > oo. The integral 
in (|3.5D is well defined in the interval [i?, oo), if f„ > 0, x > 0. Because we have chosen the 
range of z as the interval [i?, oo), we should also discuss the boundary conditions at z = R. 
We might want to follow the discussions in Case (I), that is, we require that the normalization 

■ when z ^ oo, the normalization 



conditions 
conditions 



[An - (Am) 



51) are satisfied. Because Fn{z), Gn{z 
require that 

. R' 

J I 

R 

d 

dz 



F —F* 

dz 



^ m 1 ^ n 



+ 



(A„ + \*„^){F*fin 



G„iFn) 



\R 



dz 
0. 



d_ 

^ dz 



m 1 



\R 



(3.20) 



However, it is difficult to require the solutions (|A.l| )- (|A.3D to satisfy the conditions ( 3.20]) . From 



(|3.18D , we know that are determined by the parameters ^s, ? and n. So the conditions 



( 3.20 ) impose restrictions on the parameters ^s, | and the boundary parameter R instead 



of A„. The naive numerating of parameters may mean that we can have 3 eigenstates to be 
orthogonal, because we have 3 parameters ^s, | and R. Nevertheless, such choices are difficult 
to be implemented and it seems less natural. It is more natural to regard the parameters ^s, | 
and R as the input parameters, or they should be determined by unknown physics that we 
do not consider here. In the following discussions, we will not impose boundary conditions 
further at z = Rto determine the solutions, but we will simply give these parameters by hand 
to determine the solutions. In such a choice, the normalization conditions ( |3.5| ) are not satisfied. 
So we should change to the Case (II), that is, K is matrix valued, and we try to diagonalize 
this matrix to get the action ( 3.14| ). Before doing that, we give numerical examples to show 



that only three eigenvalues are left and they are of the same order. Let | = 9/4, —s = 4. From 



Eq. ( 3.19 ), we know that only n = 1,2,3 are permitted. The eigenvalues are given by 

Ai = 3.187 UJ {n = 1), A2 = 3.833 uj (n = 2), A3 = 3.9995 lo {n = 3), (3.21) 

which are of the same order. These massive modes together with the zero mode make K and 
M to be 4 X 4 matrices. In order to get the action ( 3.14| ), K must be positive-definite. In 
Appendix we give numerical examples to show that K is positive-definite and the modified 



eigenvalues A„ in (|3.13D are still of the same order. A„ are given by 

Ao = -2.69625 u, Ai = 4.00816 uj, A2 = 4.44389 uj, A3 = 5.26792 uj. 



(3.22) 
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Here we give some interpretations for our choices of the parameters. ^ = 9/4, ^s = 4 are 
chosen to ensure that only 3 massive modes are permitted. We choose xq = = 30 as the 
boundary value in order to ensure that the modified A^, are still of the same order. We found 
that small xq, for example, xq = 1, makes different A„ to have big difference. Now we can 
understand the role of the parameters s and b in the metric (|3.15| ) . In order that we can trust 
our analysis, the condition ^ < 1 should be satisfied. From the above, we see that s appears 
in the combination ^s. Then we can always keep ^ < 1 by adjusting the value of s despite of 
the input value of ^s. While h appears in the combination so it is closely related to the 
boundary value of z. The role of a is less obvious because it appears in a more complex way. 

From the above, we notice two obvious changes: (1) The massive modes are modified to 
different size, but they are still of the same order; (2) The zero mode mixes with the massive 
modes. By this mixing, the zero mode gets mass of the same order with the massive modes. A 
strange point is that the zero mode gets a negative mass. However, it does not form problems 
in the models like [^, where only the size of the mass is relevant. In models where the sign 
of mass is relevant, we must reconsider whether it produces problems for our model. This 
new feature can supply a possibility to bypass the zero mode problem that we introduced in 
Sec. ^. In the above, we get 4 massive modes from the previous 3 massive modes and 1 zero 
mode. They can produce 4 fermion generations in 4-dimension. This is not realistic. The 
above numerical example suggests us to start with 2 massive modes and 1 zero mode. If the 
zero mode gets mass of the same order with the massive modes through mixing as the above 
numerical example, we can get only 3 generations. 

Before giving an example about this situation, we should discuss another mass source 
about the zero mode, that is, the zero mode can also get mass through coupling with a Higgs 
field on the brane sited aX z = R. Here it is appropriate to introduce the brane coupling. The 
Wilson line phase in is not well defined because the range of z is noncompact. We do not 
suggest a concrete form for this coupling. For example, it can arise from the coupling used in 
Q. Here we accept the result that the zero mode gets mass 

Ao = e w. (3.23) 

In the following discussions, we will consider two cases: (a) e ^ 0, so it is negligible; (b) 
e ~ 1, so it is comparable with the massive modes. We will give numerical examples about 
these two cases respectively. According to the same spirit with the above example, we let 
I = 8/3, = 3 to ensure that only 2 massive modes are permitted. These massive modes 
are given by 

Ai = 2.548 to {n= 1), A2 = 2.969 oj {n = 2). (3.24) 

We still choose xq = ^-5- = 30, then by the same procedure with that in Appendix we can 
get the modified A^ as 

Ao = -1.48838 cj, Ai = 3.16287 w, A2 = 3.85202 w, for e = 0; (3.25) 
Ao = 1.71708 cj, Ai = 2.59758 w, A2 = 3.19837 w, for e = 2. (3.26) 

Here we only give the results. The details are similar to that in Appendix The zero mode 
gets mass of the same order with the massive modes in both cases. The difference is that the 
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sign of the zero mode mass is opposite in these two cases. In both cases, with the help of the 
zero mode mixing with massive modes or the zero mode coupling with Higgs field, we may 
suggest a possibility to bypass the zero mode problem in Sec. In each case, we may get just 
3 generations. In models where the sign of mass is relevant, we can check which case may be 
realistic. 

In the above example, the orthogonal conditions ( |3.5D are not satisfied, because we choose 
a special metric and a special range for z. This choice induces mixing between different 
modes. In the following, we try to construct an example, which can ensure that the orthogonal 
conditions ( |3.5D are satisfied. This example can be constructed by changing the metric ( 3.15| ) 
to the following form, 



s, a,b,uj > 0. 



(3.27) 



Because a,6 > 0, this metric develops singularity^ at the point z = This may make this 



metric unrealistic. However, we find that it can satisfy the orthogonal conditions (3.5) just 
because it has such special structure. Here we let aside the problem of singularity, and focus on 
how it can satisfy the orthogonal conditions. The solutions can be given by hypergeometrical 
functions yet. 



F{z) = Cie"'''^^(e'^^ - 6)^hypergeom ( p - fj, - i^, p - p. + - 2^, — ] (3.28) 



+ C2e^'^'^(e'^'^ - 6)''hypergeom [ p + p - i', p + p + + 2p, — 



where p, p and u keep the same form with that in Eqs. ( 3.16| ). By use of the property of the 



hypergeometrical function, we see that if Re(l — 2/>) < 0, a singularity happens at z = R = 
Then if we choose the range of z to be (— oo, R], like the above example, the boundary 
conditions aX z = R = impose the conditions like (3.18). We should also make the 
solutions well behaved when z — oo. These two requirements can be satisfied by the following 
condition 



p + p — V = — n, n = 0, 1, 2, ■ 
For f < 1, Eq. ( pOj) has solut ions and n is limited by 



(3.29) 



< n < 



1 



-V 

h) 



1 



m 



LO 



(3.30) 



We give these solutions in Appendix ^ explicitly. We suppose that = 10, f = ^, in order 



that only 3 massive modes are permitted. Besides these massive modes, Eq. (3.29) also has a 
zero mode solution. These solutions are given by 



Ao = 0, Ai = 3.80 w, A2 = 4.642 w, A3 = 4.953 w. 



(3.31) 



^See Appendix ^ 
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Prom Appendix |D[ we know that when z — oo, then a; — > 0, Fn{z), Gn{z) — > 0; while at 
z = R = X = 1, then Fn{z), Gn{z) = 0. So in this example, the boundary conditions 
( |3.8| ) are satisfied. Then the orthogonal conditions ( |3.5| ) are ensured. There are no mixing 
among different modes. So in this example, the zero mode can only become massive through 
coupling with Higgs field. If the zero mode can get mass comparable to the massive modes, 
we can adjust the parameters and | to make that only 2 massive modes are permitted. 
Then we can get just 3 generations. If this zero mode gets small mass, a generation heavier 
than the SM generations is produced. We should check whether it is allowed by experiments. 
If it is objected by experiments, it will make a problem for our model. 

In Appendix ^, we suggest another example, in which the orthogonal conditions ( |3.5| ) are 
satisfied. Like the metric ( |3.27 ), there exists a singularity in the range of z we considered. 
In order to avoid the singularity, we may change the range of z. For example, for the metric 
( 3.27 ), we can choose the range of z to be where zi > In this new range, the 

singularity of the metric ( |3.27| ) at is avoided. But the orthogonal conditions ( |3.5[ ) will be 
not satisfied. We have not found a metric which satisfies the requirements: 1) it can produce 
finite generations; 2) it ensures that the orthogonal conditions ( |3.5[ ) are satisfied, and is well 
behaved in the range of z. 



4. Further Discussions 



In this section, we compare the approach adopted in |^ |5| with our setup in Sec. ^ at first. In 
the approach adopted in |5| , the authors reduced the 6-dimension spacetime to 4-dimension 
spacetime directly. Distinct from this approach, we reduce two layer warped 6-dimension 
(4-|-l)-|-l spacetime to the 5-dimension 4+1 spacetime which is still warped at the first step, 
then we reduce the 5-dimension spacetime to the physical 4-dimension spacetime. Their dif- 
ferences are the different ways that one treats the zero mode and massive modes. 

(1) For the zero mode: In the approach adopted in ^], the authors reduced the 6- 
dimension spacetime to 4-dimensions directly, and got zero modes in 4-dimension spacetime, 
so these zero modes correspond to the standard model (SM) generations. These zero modes 
get mass through coupling to the Higgs field. In the present approach, we reduce the 6- 
dimension spacetime to 5-dimensions at first, so we get a zero mode in 5-dimension spacetime. 
When one further reduces the 5-dimension spacetime to the physical 4-dimension spacetime, 
the zero mode in 5-dimension spacetime can produce a very heavy fermion in 4-dimensions, 
as illuminated obviously in [^]. It is very heavy, hence it does not correspond to the SM 
generation. However, as we discussed in Sec. |^, the zero mode can also become massive, then 
it can produce SM generation if it can get large mass. But if it gets small mass, then it 
produces a new generation objected by experimental data. This can cause a problem for our 
model. 

(2) For the massive modes: In the approach adopted in [^, ^ , the authors got the massive 
modes in 4-dimension spacetime. The massive modes are heavy Kluza-Klein (KK) particles. 
They do not correspond to the SM fermions. However, in the present approach, we reduce the 
6-dimension spacetime to 5-dimensions at first, so we get the massive modes in 5-dimension 
spacetime again. These massive modes are KK states in 5-dimension spacetime. When one 
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reduces further the 5-dimension spacetime to the physical 4-dimension spacetime, these massive 
modes in 5-dimension spacetime can produce massive fermions in 4-dimension spacetime, which 
are light and correspond to the SM fermions, as it is obvious in |^. 

These differences provide a new chance that we can address some extra issues in the present 
approach: we give an explanation for the origin of the same order bulk mass parameters, and 
give an answer for the fermion generation puzzle in the meantime. Note that it is the special 
metric Ansatz ( p. ID to supply such an explanation. The two layer structure of the metric 
enables one to reduce the 6-dimension spacetime to 5-dimension spacetime. The metric Ansatz 
for A{y) can be the AdS metric, 

A(y) = (4.1) 



By this choice, the five dimensional Dirac equation ( |2.1[l| ) is similar to that analyzed in models 
i- 

Besides, we want to address another two issues: 



(1) Whether the metric (2.1) can be the background solutions of Einstein equations? The 
metric Ansatz ( |2.1| ) has been analyzed in [|^] in high derivative gravity with matter sources and 
in |l^] with a negative bulk cosmological constant. Their solutions are not the metric which 
we suggested in Sec. |^. Here we consider a minimum coupled scalar-gravity system in order to 
investigate whether the metric in Sec. ^ can be realized. This is a simple and convenient way. 
The action is given by 



d'xdydzy/^ilM^R] + / d^xdydz^/^\ ^g^'^^V M(py + Vi(p) } , (4.2) 



S 



in which V{(f)) is the potential term for scalar field. Supposing that the metric Ansatz 
and that (p only depends on z, we get the following equations, 



1 



yy 



4 a2 



IQB^'^BI + 4A~^A 



yy 



B- 



+ AB-'B, 



+ 2A-^Al 
dV{(j)) 



4M4 
1 

4M4 
1 

4M4 
0, 



-B' 



-B' 



'-B- 
2 

2 



^<t>l + v{<t>) 



B^i\B 



-2 i,2 



(4.3) 
(4.4) 
(4.5) 
(4.6) 



in which A,, 



dA 
dy ' 



B. 



dB 
dz ' 



dz ' 



Obviously, A(y) should be of the form A(y) 



ky+c 



m 



which k, c are constant. A minimum coupled scalar-gravity coupled system has been analyzed 
in Ijl^ in five dimensions. The result is that there always exists appropriate form of V{4>) 
to ensure that the metric have solutions, and V{(p) and the metric can be expressed with a 
superpotential. The similar result applies to the above system, that is, for any B{z), there 
exists appropriate V{<p), which makes Eqs. ( [4.3D -(^) satisfied. However, the present system 
is more complex, and it is difficult to express the solutions with a superpotential. Of course, 
in order to get the solution of the metric adopted in Sec. ^, the boundary conditions must be 
adopted appropriately. 
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(2) Whether the examples used in Sec. ^ are acceptable physically? The metric (|3.15 ), 
( |3.27 ) and ( [E.l| ) are similar to that analyzed in |19|, that is, they are both asymptotically flat. 
They are also both noncompact, and both have infinite volume |2^]. Especially for the metric 
( |3.27 ) and the metric ( [E.l| ), there is singularity |21] in the range of z that we choose. So it needs 
further work to investigate whether they are acceptable physically. The similar problem exists 
for the metric ( p.22| ), in which the sixth space dimension is timelike. As being emphasized 



in [22 1, the violations of casuality and probability give stringent restrictions on the timelike 
dimension. It also needs further work to investigate whether it is acceptable physically. We 
have not found a metric, which is finite volumed, as in the models and can produce finite 
generations simultaneously. 



5. Conclusions 

Now we summarize the main points in our work. In this paper, we try to explain the ori- 
gin of the same order bulk mass parameters, and give answers to the generation replication 
puzzle simultaneously. The fermion masses are of hierarchy structure in 4-dimension space- 
time. It seems that it is difficult to interpret them as the eigenvalues of a Schrodinger-like 
equation. However, the hierarchy structure can be reproduced with the bulk mass parameters 
in 5-dimension spacetime. The 5-dimension mass parameters are in the same order, as have 
been shown in many papers |^, ^J. This interesting feature supplies a chance to interpret 
the 5-dimension mass parameters, which are of the same order, as the the eigenvalues of a 
Schrodinger-like equation. Supposing that the six dimension spacetime metric has special two 
layer (4+l)+l structure, we can reduce the 6-dimension spacetime to 5-dimension spacetime at 
the first step. We find that the bulk mass parameters are the eigenvalues of a Schrodinger-like 
equation. Hence the same order mass parameters emerge naturally. However, the problem is 
that the number of eigenvalues is infinite generally, which leads to infinite light generations. We 
suggest several approaches to deal with this problem. Obviously, this problem arises from the 
fact that in the conventional Kluza-Klein (KK) decomposition, one gets infinite KK particles 
generally. However, as in the example given by Madore (2^, in the noncommutative geomet- 
rical background, and by the choice of the internal structure, the modification of KK theory 
gives rise to finite spectrum of particles. Therefore, it is possible to overcome the difficulties 
in our work by use of the noncommutative geometry. It will require modifying the framework 
in Sec. |. We hope we can address these issues in the future. 
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A. Explicit solutions for massive modes and zero mode: metric I 

In this appendix, we give the solutions of Eq. ( 3.16| ) for massive modes under the conditions 
(3.18) explicitly. The solutions should be well behaved in the range oo). Let x = the 
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solutions are given by 



N2 

UJ 



ai X 
7i X + 1 



^ _ 2a3 X _^ 03(03 + 1) 



73 a; + 1 73(73 + 1) \x + l 



X 



(A.l) 
(A.2) 

(A.3) 



in which 



7n = 1 - 2/i„, n = 1,2,3, •• • 

m / a\ 



CO / \ UJ 



m \2 /a\2 / A„ 



(A.4) 
(A.5) 

(A.6) 
(A.7) 



We have dropped the hyper geometrical function after the coefficient C2 in ( ^.161) , because it 
is divergent under the conditions ( p.l8 ) when z — > 00. The solutions for Gn{z) can be gotten 
from ( |2.14| ) as 



Gn{x) 



ijj 
An, 



m x+j- , ^ d 

— S—^Fn{x) - X — Fn{x) 

UJ x + \ ax 



The zero mode solution is given by 

Fo{x) = 0, Go(x) 
When z — > 00, X = > 00 



/]Vo 



^ x"'^", Gn{x) (— S + l/„) X 

Fo(x) = 0, Go(x) ^x-t^ 



(A.8) 

(A.9) 

(A.IO) 
(A.ll) 



They are all well behaved when f„ > 0. In terms of x, the integral in Eq. ( l3.2| ) can be rewritten 
as 

f 1 f dx 

Kmn = / dz {F:,Fn + Gl,Gn) = - / — {KFn + G^^Gn) ■ (A.12) 

J J X 

They are also well behaved in the range [R, 00) when f„ > and x > 0. In the numerical 
examples we give in Sec. 0, the conditions f „ > and x > are always satisfied. 
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B. Numerical examples for finite generations 



In order to get numerical results, we need to input the parameters ^s, f and R. From the 
solutions in Appendix A, we know that it is enough to input value for xq = 

Let '^s = 4, ^ = 9/4 as in Sec. ^, then for massive modes, only n = 1, 2, 3 are permitted. 
We further designate xq = = 30. We normalize the solutions for massive modes and zero 
mode according to (|3.5| ) as 



J dz {F*Fn + GlGn) = 1, n = 0, 1, 2, 3. 



(B.l) 



These conditions determine the normalization constants. Then the matrix K is determined to 
be 



K 



( 1 0.8797 0.7144 0.3503 \ 
0.8797 1 0.9399 0.5173 
0.7144 0.9399 1 0.6433 

\ 0.3503 0.5173 0.6433 1 / 



(B.2) 



The indices of K are determined according to ( ^.2| ) as 



Kmn = y (^m-^n + G*^Gn) , m, n = 0, 1, 2, 3. 



(B.3) 



K can be diagonalized as 



K = V^KV, A = diag(3.0656, 0.707478, 0.211924, 0.0149942), 

/ -0.491317 -0.555244 -0.544655 -0.391998 \ 
-0.488159 -0.224439 0.0638874 0.84098 
-0.651525 0.221179 0.626183 -0.366729 
\ 0.309555 -0.769683 0.554224 -0.0678292 / 



V 



(B.4) 



where V'^ means the transpose of V . This example shows that K is positive-definite, as we 
expected. 

A„ is given by ( 3.21| ) in Sec. ^. M is determined by ( ^I^ ) as 



dz 



[Prn^n + G*^Gn) 



An + 



, m, n = 0, 1, 2, 3. 



(B.5) 



Given the parameters, M is determined to be 



/ 1.402 1.3693 0.7006 \ 



M 



1.402 



/1463 
12 



3.2993 1.8589 



1.3693 3.2993 2.5196 
V 0.7006 1.8589 2.5196 / 



(B.6) 
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Following the procedure in Sec. 0, we can get M as 



M 



U 



A 



( 2.73442 -1.35109 -2.28483 2.23697 \ 
-1.35109 3.00477 -1.47008 1.60816 
-2.28483 -1.47008 2.19631 1.7851 
\ 2.23697 1.60816 1.7851 3.08822 / 

/ 0.671378 0.216423 -0.209452 0.677158 \ 
-0.442926 0.350635 0.637256 0.524191 
0.256197 -0.814717 0.495127 0.159525 
\ -0.536126 -0.407986 -0.552162 0.491155 / 
diag(5.26792, 4.44389, 4.00816, -2.69625). 



(B.7) 



C. Property of hypergeometrical function 



We cite a theorem about the hypergeometrical function -F(a, /3; 7,^). 

Theorem 2.1.2 The series F(a^ (3; 7, with |^| = 1 converges absolutely if Re(7— a— /?) > 
0. The series converges conditionally if ^ = e*^ 7^ 1 and > Re(7 — a — P) > —1 and the series 
diverges if Re(7 — a — P) < —1. 



D. Solutions for zero mode and massive modes: metric II 



For the metric ( |3.27| ), the solutions for massive modes well behaved in the range {—oo,R] are 
given by 



Fi{x) 
F2{x) 



x^'^{l-x)f 
xf'^l-xy 



-X 



71 
2p2 



, P2{P2 + l) 2 

X H 7 —X 



72 72(72 + 1) 



(D.l) 
(D.2) 



in which 



l + 2nn, n = l,2,3,' 



In 



Pn = P + Hn + I'n, 



An 
UJ 



(D.3) 
(D.4) 

(D.5) 
(D.6) 



In the above expressions, we have defined x = p is defined as in Sec. ^ are determined 
by (3.29), and n is limited by ( p. 30 ). We have dropped the hypergeometrical function after 
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the coefficient Ci in ( p.28 ), because it is divergent when z 
are determined by 



-oo. The solutions for Gnix) 



Gn{x) = ^ 



— S ^-Fn{x) - X—Fnix) 

u) X — 1 ax 



(D.7) 



The solutions ( D.1 )-( |D?7|) are well behaved when ;U„ > 0, p > 1. They are satisfied in our 
numerical example in Sec. |3[ 

We find that Eq. ( ^.291) also has a well behaved zero mode solution. This zero mode 
solution is given by 



Fo{x) 



x'^'il-xY, Go(x) =0. 



(D. 



It is consistent with what we get from ( 2.16| ). 



E. Another metric example for orthogonality 



In this appendix, we suggest another metric which ensures that the the orthogonal conditions 
( |3.5| ) are satisfied. This metric is given by 



ujz - a 

B[z) = s -, s,a,b,uj > L). 

ojz — b 



(E.l) 



We suppose a,b > and a > b here. The conditions a, & < and a < b work well also. But 
the conditions a < 0,6 > do not work. Eqs. ( |2.17 ) and ( p.l8| ) can be solved by confluent 
hypergeometrical function (or Kummer's function). When 7 = 2(1 + z^), 1/ = ^s{a — b) is not 
integer, the solutions are given by 

F{z) = Cie-'^('^^-*)(u;z - 6)i+^hypergeom (a; 7, 2/i(wz - 6)) 

+ C2e-'^('^^-*)(wz - 6)i+'^+i-Thypergeom (a + 1 - 7; 2 - 7, 2/x(wz - b)) , (E.2) 



in which p = y (^-s)^ — (^)^ and a = 1 + u — ^s^. We choose the range of z to be [R, 00), 
where R = ^. The confluent hypergeometrical function F(a; 7, ^) ~ when ^ — > 00. In order 
to make the solutions well behaved when z ^ 00, the confluent hypergeometrical function 
must be cut off to be a polynomial by the requirement 



a = 1 + 
Then A„ are determined to be 



m u 

— s— = —[n—l), n=l,2,---. 



-sll 



{u + ny 



to. 



(E.3) 



(E.4) 



By ([E.4|), we know that 



m (2z^ + l 

— T 



< A < ( — s ) oj. 



(E.5) 



- 19 - 



We should also require that the solutions are well behaved at 2; = i? = ^. This condition 
requires the solutions further to be 



Fniz) = ^^e-f'^''-^\x - 6)i+'^hypergeom(-n;7,2/i(x - b)) . (E.6) 
v^i 



Here we define x = coz, and Gn{z) are determined by 



Gn{x) = ^ 
An 



m X — a d 

— S rPnix) - —Fn{x) 

00 X — ax 



There is also a well behaved zero mode solution, which is given by 



(E.7) 



Fo(x) = 0, Go(x) = ^e-^^-(x - by. (E.8) 

We have supposed a > b. So when z — > 00, x — > 00, Fn{z), Gn{z) 0; while Fn{z), Gn{z) = 
aX, z = R, X = b \i V > 2. So the orthogonal conditions ( |3.8D can be satisfied. 

According to the analysis in 1^], if we let Aq = ms to be the lightest generation of SM, 
there exist infinite heavier generations corresponding to A < ms, in which A = is the heaviest 
generation. However, there exists a problem in this case: the lighter generations approximate 
to be continuous, which conflicts with the experimental fact. Therefore, special boundary 
conditions must be adopted to remove the reductant generations. We hope that there are only 
finite generations left with the help of the special boundary conditions. But it seems that it is 
difficult to impose such boundaries naturally. 



F. Ricci scalar curvature for metric 



Ricci scalar curvature for metric Ansatz ( |2.1| ) is given by 



R = - [10(5-3^,, + B-^Bl) 



+ (8A- 



^Ayy + 4A 



-4 a2 



(F.l) 



in which A,, 



dA D 
dy ' 



^ . From the second term in 



, we know that there is singularity 
^. At the point z = ^ and z = 
the metric ( J3.27| ) and (E.l) are not well defined respectively, but the Ricci scalar for them are 
well defined. So it needs further work to determine whether they are true singularities. 



for the metric ( |3.27] ) at z = and the metric (E.l) at 
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